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Abstract 

Gravity is treated as a stochastic phenomenon based on fluctuations of the 
metric tensor of general relativity. By using a (3+1) slicing of spacetime, 
a Langevin equation for the dynamical conjugate momentum and a Fokker- 
Planck equation for its probability distribution are derived. The Raychaud- 
huri equation for a congruence of timelike or null geodesies leads to a stochastic 
differential equation for the expansion parameter in terms of the proper time 
s. For sufficiently strong metric fluctuations, it is shown that caustic singular- 
ities in spacetime can be avoided for converging geodesies. The formalism is 
applied to the gravitational collapse of a star and the Friedmann-Robertson- 
Walker cosmological model. It is found that owing to the stochastic behavior 
of the geometry, the singularity in gravitational collapse and the big-bang 
have a zero probability of occurring. Moreover, as a star collapses the prob- 
ability of a distant observer seeing an infinite red shift at the Schwarzschild 
radius of the star is zero. Therefore, there is a vanishing probability of a 
Schwarzschild black hole event horizon forming during gravitational collapse. 
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I. INTRODUCTION 



Attempts to formulate a quantum theory of gravitation have not met with success to 
date The main attempts have been based on: 

1. Perturbation theory based on expansions about Minkowski spacetime; 

2. Canonical Hamiltonian formalism; 

3. Path integral formalism; 

4. String theory. 

The first approach became untenable for two reasons: (a) The series expansions were 
not renormalizable for gravitational couplings to matter and beyond the first order for vac- 
uum gravity, (b) The method chose a special frame of reference, thus failing to meet the 
fundamental requirements of general relativity. Approach (2) ran into problems with the 
Hamiltonian constraints, since one set of constraints was non-polynomial and could not be 
solved to yield the true gravitational degrees of freedom. Attempts to resolve this problem 
have met with some success, but fundamental difficulties still exist that prevent a satisfac- 
tory formulation of quantum gravity Q]. In particular, it is difficult to meet the necessary 
criterion that the quantum gravity theory must have a well-defined classical general rela- 
tivity (GR) limit ||^. Approach (3) has met with the unresolved difficulties of defining a 
meaningful measure in the spacetime manifold, and the need to use an Euclidean spacetime 
to define the path integral. In general, it is not possible to return to curved spacetime with 
a local Lorentzian signature by a Wick rotation, and guarantee that the classical GR limit 
exists. Finally, string theory has not succeeded in describing an easily interpretable so- 
lution to quantum gravity, for the theory must be formulated as a non-perturbative scheme 
and potential problems with causality and non-locality prevent the use of standard canonical 
Hamiltonian methods and other field theory techniques based on perturbation theory. 
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Sakharov and Jacobson have formulated effective gravitational theories. Sakharov 
viewed spacetime as an elastic bulk system with an elasticity constant. However, this theory 
had difficulty yielding Newton's and Einstein's gravitational theories in a well-defined limit. 
Jacobson pictured Einstein's field equation as an equation of state and the gravitational 
variables were thermodynamic variables, so that the Einstein equation was derived from 
the proportionality of entropy and the black hole horizon area together with the relation 
5Q = TdS connecting heat Q, entropy S", and temperature T. In this approach, it is not 
considered appropriate to quantize the classical gravitational fields as is done for the case of 
the electromagnetic fields. 

It is well known that an analogy between classical statistical mechanics and Euclidean 
quantum field theory exists. This is based on the fact that the path integral measure 
D0exp(— S'£;[0]) can be identified with the Boltzmann probability distribution of a statistical 
system, where Se[<P\ denotes the field theory action, which has been Wick rotated to a 4- 
dimensional Euclidean space. The Euclidean Green functions are the correlation functions 
of a statistical system in thermal equilibrium. A stochastic quantization of field theory has 
been formulated by Parisi and Wu [0,^. A new 'fictitious' time coordinate was introduced 
with respect to which all fields evolve, namely, the field was supplemented with the 
extra coordinate t, such that (f){x,t), where ) and the fictitious 

time t should not be confused with the Euclidean time x". The evolution of (f){x,t) is 
determined by the Langevin equation: 



where C,t{x, t) is a Markov process (Gaussian white noise). The basic assumption of stochastic 
quantization is that in the limit t ^ oo equilibrium is reached, and that the equal time 
correlation functions ^[r/] become identical to the corresponding quantum Green functions: 



where the right hand-side denotes the standard vacuum expectation value in Euclidean 




(1) 




(2) 



space. 



3 



The stochastic quantization method has been used as a way to quantize the gravitational 
field by Rumpf 0. However, since to any ordinary Feynman diagram there corresponds a 
number of stochastic diagrams with the same topology, and it can be shown that in the 
equilibrium limit these stochastic diagrams exactly sum up to the corresponding Feynman 
diagram, then the standard problems of perturbative quantum gravity will still prevail. 
Moreover, there is no indication that the method alleviates any of the problems of non- 
perturbative quantum gravity. 

Over the past few years, theoretical and experimental studies have shown that a non- 
linear system in a "noisy" environment displays surprising behavior that does not conform 
to common intuitive experiences. Self-organization in non-equilibrium systems coupled to 
fiuctuating environments exhibits a much richer behavior than is possible in a determinis- 
tic situation. Transition phenomena show features similar to equilibrium phase transitions 
and to transition phenomena encountered in non-equilibrium systems influenced by external 
constraints, e.g., the Benard instability and the laser transition. Phase transition phenom- 
ena and chaotic behavior can be induced by environmental randomness |]TD|]. Statistical 
mechanics can display two kinds of systems: (1) a fine-tuning of parameters that leads to 
critical points and phase transitions and (2), self-organization which can often describe crit- 



ical systems with structure spread out over every available scale [|rT|. These properties can 



be applied to an understanding of our universe at both the large and small scale |12 



In quantum gravity there is only one scale that contains Planck's constant h, namely, 
the Planck length, Ip = ^/fiG, so that quantum gravity describes a strongly coupled phase 
in which there are no correlations on large scales. However, space and time can themselves 
be consequences of critical phenomena in which there occur fiuctuations with correlations 
on a scale much larger than the Planck scale. Thus, there can exist a gravitational system 
which is weakly coupled and critical showing self-organization behavior. Such gravitational 
phenomena can help to remove fine-tuning behavior encountered in cosmology e.g., those 
associated with the initial conditions of the universe and the cosmological constant. 
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It is plausible that a nonlinear gravitational system coupled to a stochastic, fluctuating 
spacetime geometry can lead to significant changes in the macroscopic behavior of the sys- 
tem, for the fiuctuations could be of order and not where V denotes the characteristic 
volume of the system. Another way to ask this question is: do nonlinear gravitational sys- 
tems, coupled to a rapidly fiuctuating spacetime geometry, always adjust their macroscopic 
behavior to the average properties of the spacetime geometry or do there exist situations in 
which the system responds to the randomness of the spacetime geometry, displaying behavior 
forbidden under deterministic conditions? 

Attempts have been made to treat the geometry of spacetime in a semi-classical way by 
writing Einstein's gravitational field equations in the form: 

G,, = 87:G{T^,), (3) 

where (T^u) denotes the expectation value of the stress-energy tensor operator T^u in a given 
quantum state. In order to reconcile the left hand-side of (0) with the right hand-side, it 
is necessary to make the gravitational field change in a discontinuous, acausal manner, a 
behavior unacceptable from the classical point of view. 

We shall develop a probabilistic interpretation of spacetime by defining the metric tensor 
g^u as a stochastic random variable, which is defined by a set of possible values and a prob- 
ability distribution P{g) over this set, with P{g) > and normalized to unity on its range, 
where g denotes the metric tensor g^,^ and P{g) is a scalar obtained from the contraction 
of a tensor quantity invariant under diffeomorphism transformations. The probability that 
g takes a value between g and dg is given by P{g)dg. We shall use the (3+1) foliation of 



spacetime in GR |[13 Hlq] and the canonical formalism to derive a Langevin equation for the 
conjugate canonical momentum, which is treated as a random variable. A Fokker-Planck 
equation for the time evolution of this stochastic variable and an equation for the probability 
density are derived. In this formalism, we treat both the gravitational field and the matter 
in a probabilistic fashion, so that Eq. (|^) now becomes 

(G^,) = 8nG{T^,), (4) 



and there is no logical conflict between the two sides of the equation, although the Gf^u and 
the Tf^i, are not treated as quantum mechanical operators but instead are considered to be 
functions of random, stochastic variables. 

In Sections 11 through V, we introduce basic notation for stochastic theory and probabil- 
ity as applied to GR, and in Section VI, we derive the Langevin equation for the conjugate 
momentum variable in (3 + 1) foliated spacetime and a Fokker-Planck equation for the prob- 
ability density. In Section Vll, we use the Raychaudhuri equation to derive a stochastic 
differential equation for the expansion parameter 6, and investigate the stochastic behavior 
of a converging congruence of geodesies. In sections Vlll and IX, we apply the method to 
gravitational collapse and cosmology and we conclude, in Section X, with a summary of the 
results. 



II. METRIC FLUCTUATIONS AND STOCHASTIC PROCESSES 



Wheeler |]T^ suggested that quantum fluctuations in the metric tensor of spacetime, Qf^i,, 



should occur at the order of the Planck length: Ip = (hCy^'^ = 1.6 x 10 ^"^cm. These 
fluctuations are of order: 

Ag ~ ^, (5) 

where L is the spatial dimension of the region in a local Lorentz frame of reference. Fluctu- 
ations in the curvature tensor in space are of the order: 

For a domain of order 1 cm, the curvature fluctuations are 



AR ~ 10~^^cm-^ (7) 

and for subatomic domains the metric fluctuations are Ag ~ 10^^°. These are so small that 
we can ignore them for experimentally accessible physical observations. 
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Let us postulate that the spacetime geometry is a physical system. In classical GR, it 
is assumed that the spacetime manifold is smooth down to zero length scales. Since 
all known physical systems possess "noise" at some length scale, the former assumption 
would seem to be unrealistic. We shall also postulate that the spacetime geometry has 
an (unknown) underlying microscopic structure and that the subsystems of this structure 
undergo stochastic fluctuations at a certain length scale, which can produce macroscopic 
correlation lengths due to cooperative, self-organized behavior when the subsystems couple 
to matter. 

If we treat the geometry of spacetime as a fluctuating environment, then we can consider 
the gravitational constant G as a control parameter with the decomposition: 

Gt = G + a^t, (8) 

where G is the Newtonian gravitational constant ( G = 6.673 x lO^^g^^ cm^ s^^ ), and is 
Gaussian white noise with E{C,t} = where E denotes the expectation value; a measures 
the intensity of the geometrical fluctuations of the metric. Thus, G will have a bell-shaped 
curve distribution, peaked at the average value of G. The central limit theorem assures that 
the fluctuations in G are Gaussian distributed. The spacetime fluctuations of the metric are 
of order: 

^9 ~ — (9) 
Let us describe the gravitational system by a set of state variables X obeying the equation 



dtXir,t) = fGiX{r,t)), (10) 

which describes the local evolution of X in space r and time t. The state variables X{t) 
are constructed from the metric tensor g, its derivatives and possible curvature quantities, 
and the functionals fciX) depend on the control parameter G. When the symmetry of the 
gravitating system is broken, self-organization sets in and the disorganized state becomes 
unstable. We assume that there exists a time independent solution X satisfying 



foiX) = 0, (11) 
corresponding to an unorganized state. We expand X about X in small perturbations: 

X{r,t)^X + x{r,t), (12) 
where \x/X\ <^ 1, and the perturbations x{r,t) satisfy the equations: 

dtXi^^BikXk, (13) 

k 

with the solutions: 

^^"(0 exp(a;fci). (14) 

k 

The hfe-time of gravitational disturbances is of order 

^ (15) 



' macro t-> / \ 

Re(a;k) 

where Tmacro is the macroscopic time scale of evolution of the system. The Re(a;fe) must be 
negative for a gravitational state X which is asymptotically stable. Gravitational transitions 
can occur at the point G = GcSit which Re(a;/j) changes sign from negative to positive values, 
and a new solution can occur at this point of bifurcation associated with a gravitational order 
parameter which vanishes at Gc- 

The correlation time Tcor is a measure of the rapidity of the random fluctuations of the 
spacetime geometry. For a stationary gravitational process it is defined by 

Tcor = / C{T)dT, (16) 



C(0) Jq 

where the normalized correlation function is given by 

~ E{5gt5gt+r} , . 

""^^^ ^ E{Sgn • ^''^ 

For a typical exponentially decreasing correlation function, 

C{T)^C{0)eM-ir), (18) 



where 7 is a constant, Tcor = 1/7. Spacetime will be characterized by rapid fluctuations with 

^cor ^macro- ('^^) 

For a given 3-geometry ^^''Q, including geometrical noise in the gravitational equations 
leads to the formal Langevin |TP| or stochastic differential equation (SDE): 



dtXigt)=fGA9t), (20) 

where gt denotes the random metric tensor. Assuming for the moment that the external 
parameter Gt enters in a linear way, then ( ^Of ) takes the form: 



dtX = f{gt) + Gh{gt) + aiMgt)- (21) 

The extreme irregularity of Gaussian white noise means that the time derivative of the 
metric is not well defined. However, the standard theory of stochastic processes can handle 
this difficulty by defining the equivalent integral equation: 

Xt = Xo + FcMdy + aj^ h{gy)iydy, (22) 

where 

FGA9t) = f{9t) + Gh{gt). (23) 



Two definitions of the stochastic integral / h{gy)^ydy have been given by Ito ||20[ and 



Stratonovich pT|. Both of these definitions are based on the concept that integration of 
Gaussian white noise yields Brownian motion, which we denote by Wt, i.e., C,t = dfWt. 
Then, 

/ h{gy%dy = J h{gy)dWy. (24) 
Both the Ito and the Stratonovich definitions of the integral lead to a consistent calculus. 
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III. ELEMENTS OF METRIC FLUCTUATION PROBABILITY THEORY 



For a manifold Ai in 4-space, we assume the existence of a probability triple {Q, A, P) 
consisting of a sample space fi, a field of events A and a probability measure P |^^. Only 
Vt and A are used to define a random variable in M.. The sample space f2 is the ensemble 
of elementary outcomes of spacetime measurements in M., labeled uj : uo eVL. We take it as 
given that the number of elementary outcomes in the manifold M. is uncountably infinite. 

The second element of the probability triple of is ^ (the a field) of events. The set 
A satisfies the following properties: 

(i) A contains the certain event, E A. 

(ii) It contains the empty set 0: (p E A. 

(iii) If A is an event, then so also is the complement A = — A. 

(iv) If A and B are events, so is their union and intersection. 

The metric tensor g is defined as a random variable by the requirement that it is a 
function from the sample space VL into some state space: g : Q ^ TZ. We also have 

A = X-\{-oo,g])eA, Wgen, (25) 

where 

X-\B) = {uj\g{uj) e B}, Ben. (26) 

We define the third element P as the frequency of occurrence of event A. Thus, on the 
interval [0, 1] P : A ^ [0, !]• We have by definition: 

P(0) = 0, (27a) 
P{n) = 1, (27b) 

An G A, An Hn^m Am = (f) ^ P (u^^.An) = £ P(A„). (27c) 

^ ^ n=l 

The distribution function induced by the random metric variable g is defined by 

Fg = Pi{cu\giu;)}) = Pig), g eU. (28) 
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The probability density p{g) is defined by 

F{9) = fp{y)dy, (29) 

where the integral is defined as a functional integral over the metric g with a given measure. 
In the following, it is understood that the distribution function Fg and the probability 

density p{g) transform as tensor quantities under diffeomorphism transformations of the 
coordinates. 

In the standard way, we define the expectation E{g} of the metric tensor by 

E{g}^ [ gdP= [ g{uj)dP{uj)^ f g{u;)P{du;), (30) 

and we have 

E{g} = / gp{g)dg, (31) 
Jn 

where E{g} is a tensor quantity. The rth moment of g is 

E{g^} = [ g'-{uj)dP{u) = f g^dF{g) = f g'p{g)dg. (32) 

The mean deviation of the random metric function is 

a' = E{{5gf} = E{ig - mf} = ! {g - mfdF{g), (33) 

Jn 

where m denotes the mean value of g. The number and location of the extrema of the 
probability density are 

5.P(^)L=.. = 0, (34) 

where dgp{g) denotes a functional derivative with respect to the metric g. The most probable 
states are given by the maxima: 

dggP{9)\g=g^ < 0. (35) 

We shall denote stochastic processes for the metric tensor in a 3-geometry ^^^Q by gt 
and deterministic metric processes by g. The set {gt, t & 9} of real-valued random metric 
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variables, gt : {ft, A, P) — > {TZ, B), describes a random function with the index 9 and a state 
space TZ. The index parameter will be the time and the index set 9 is either the real line TZ 
or the non-negative half-line, if the process started at t = 0. 

The standard Wiener process Wt has the intial condition: Wq — 0. The hierarchy of 
probability densities for a one-dimensional fluctuating metric is 

p{g, t) = (27ri)-^ exp(-//2t) = n{g, t), (36a) 



p{gi, ti, gm, tm) = n{gi, ti)n{g2 - gi, t2 - ti)...n{gm - gm-i, tm - tm-i)- (36b) 

Since the Wiener process is a Gaussian process, we define the finite dimensional proba- 
bility density as Gaussian: 

Pg{9) = [{2nrdetM]-'/^ exp[-l/2{g - mf M-\g - m)], (37) 

where 

g^"" = {9i,--:9n): m*"" = (mi , . . . , m„) , (38) 

and M is a positive definite n x n matrix. 

A stationary stochastic process in a 3-geometry ^^^Q has a probability density that is 
invariant against shifts of time: 

p{gu h; gn: t^) = p{gi, ti + 1; g^, tn + t). (39) 

Then, we have 

p{9,t)^Ps{9), (40) 

and 

E{9t}= / gp{g,t)dg= / gps{g)dg = m. (41) 
Jn Jn 
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Moreover, the two-dimensional probability density satisfies 

Pigi, h; g2, t2) = p{gi, 92; ^2 - h), (42) 

so it only depends on the time difference t2 — ti. The correlation function, defined by 

E{6gt^6gt^} = {gi - m){g2 - m)p{gi,ti] g2,t2)dgidg2, (43) 
Jn. Jn 

also only depends on the time difference t2 — ti- 

The Gaussian white noise is very irregular but is a useful model for rapidly fluctuating 
phenomena. Since it has no continuous sample paths and it has infinite total power: 

S = 2 (aV27r)c/z/ = 00, (44) 
Jo 

it does not occur in nature. In the idealization of 5-function correlated metric fluctuations, 
the gravitational system will be described by an SDE or Langevin equation of the form (plD, 
where a denotes the intensity of the Gaussian white noise. 

A random process for the metric fluctuations is a Markov process if when the present state 
of the process gt is known, then any additional information on its past history is irrelevant 
for the prediction of its future evolution. The only Markov processes which have continuous 
and differentiable sample paths are the deterministic ones associated with classical gravity 
theory, given in our ^^^Q by 

dtX{git)) = figit)), (45) 

with X(„ = J, where the only random element is the initial value I. For every positive e and 
every function of the metric, f{g, r), there is a drift, 

\imt-.r p{y,t\g,T)dy = f{g,T). (46) 

For every positive e, there is a function h{g, r) called diffusion, which satisfies 

limt^r {y - g fp{y,t\g,T)dy = h'^{g,T). (47) 

We shall also assume that there exists a positive 6 such that the condition: 

lim^^, f \y- g\^+'p{y, t\g, r) = 0, (48) 
Jn 

is satisfied. 
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IV. THE FOKKER-PLANCK EQUATION AND STATIONARY PROBABILITIES 



The transition probability density p{y,t\g,T) of the metric diffusion process gt satisfies 



the Fokker-Planck equation p3| , p4| (FPE) or Kolmogorov forward equation |T7]: 

dtp{y,t\g,T) = -dy[f{y,t)p{y,t\g,T)] + ^dyy[h'^{y,t)p{y,t\g,T)]. (49) 

Here, it is assumed that the functional derivatives occurring in (|^) exist and are continuous. 
As in the previous sections, we have suppressed tensor indices. For time homogeneous 
Markov metric processes, the Fokker-Planck equation is given by 

dtP{y,t\g) = ~dy[f{y)p{y,t\g)] + ^dyy[h'{y)p{y,t\g)]. (50) 

Here, the drift and the diffusion are time independent. For a one-dimensional system, the 
FPE for the probability density becomes 

dtp{g,t) = -[djig,t)p{g,t)] + ^dyy[h\g,t)p{g,t)]. (51) 

An important and useful solution of the FPEs can be obtained for stationary random 
metric processes. We expect that a gravitational system subjected to metric fluctuations 
for a sufficiently long time will settle down to a stationary behavior. This means that as 
time approaches infinity the system will reach a state for which the probability density 
Ps{g) has a shape that does not change with time, i.e., the probability to find the system 
in the neighborhood of a particular state becomes time independent. However, the sample 
paths gt{uj) will in general not approach a steady-state value gsi^j), so that the state of the 
system still fluctuates. However, these fluctuations are such that gt and gt+r have the same 
probability density, ps{g)- 

The stationary solution ps{g) of the FPE satisfies: 

dtpig, t\go, 0) + dylig, t\go, 0) = 0, (52) 

where 
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2 

/(^/,t|^7o,0) = f{g)p{gA%.^) - \dg[h\g)p{g,t\g,,Q)l (53) 

and we have for simplicity considered the time homogeneous case. The stationary FPE is 
then given by 

dgls{g) = 0, (54) 
which imphes that Is{g) = constant for g G [(71,(72]- 



The solution of the stationary FPE equation reads |T0| 



where C denotes a normalization constant. When the boundaries of the gravitational system 
are natural, i.e., regular with instantaneous reflection imposed as a boundary condition, then 
there is no flow of probability out of the state space and / = 0. By using the Ito prescription, 
we obtain 

In this case, the normalization constant C is given by 



V. NONLINEAR GRAVITATIONAL STOCHASTIC SYSTEMS 

If the gravitational equations are linear in the control parameter G, then we can apply 
directly the stochastic methods to solve for the probability density functions. However, we 
shall find that in certain applications the gravitational constant will enter the phenomeno- 
logical equations in a nonlinear way. Consider a phenomenological equation of the form: 

dtXit) = f{X{t)) + P{G)h{X{t)), (58) 
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where /3 is a nonlinear function of the gravitational constant G. For the geometrical fluc- 
tuations, we replace the constant control parameter G by a stationary stochastic process 
Gt = G + Ct and form the SDE: 

dXt = [f{Xt) + I3{G + Ct)h{Xt)]dt. (59) 

Since nonlinear operations on generalized functions, such as the Dirac (5-function, cannot 
be given a well-defined mathematical meaning, we cannot use a white-noise Gaussian ap- 
proximation by setting Gt = G + a^f However, methods exist which lead to satisfactory 



approximations to the white-noise idealization p5| , p!0 . 

We shall assume that Q is "colored" noise and that in applications the geometry of 
spacetime varies on a much faster time scale than the gravitating matter system coupled 
to it, so that (t is a process with a short correlation time. Let us define the geometrical 
fluctuations by the process rjt: 

Vt = PiG + Ct) - E{PiG + Ct)} = P{G + Ct) - m{G, a'). (60) 

Then we can write Eq.(|5SD as 

dXt = [fiXt) + m{G, a'')h{Xt)]dt + h{Xt)r]tdt. (61) 

We can "speed up" the stochastic process by writing 

Vt = Vt/e^ = f3{G + Ct/eO - m(G, a'), (62) 

where e = Tcor is a small parameter. Then we have 

dX't = [f{X't) + m{G, a'^)h{Xl)]dt + -h{Xt)r]ldt, (63) 



and 



dCt = -\Ctdt+-dWt. (64) 



e2 e 



Rescaling time, t t/Tcor-, the FPE has the form: 
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dtp{x,y,t)=^^ + ^ + Dsy{x,y,t), (65) 



where 



^2 

Di = dyV + —dyV, (66a) 

D2 = -ydM^), (66b) 

= -d.[f{.x) + m(G, a2)/i(a;)]. (66c) 

The correct SDE in the white-noise Gaussian hmit is given by 

dXt = [f{Xt) + m(G, a'')h{Xt)]dt + ~ah{Xt)dWu (67) 

where 

= 2 / p5(y)[/3(G' + y) -m(G',<7^)]^(G' + y)dy, (68) 
Jn 

and P5(x) is the stationary probabihty density, which has the form: 

^^(") = ^^"p(^/^;(§^^4 

where 

F{x) = /(x) + m(G', a'^)h{x) (70) 

and 

t;(x) = h{x)[l + egi(x) + e'g2(a;) + ...]. (71) 

VI. HAMILTONIAN FORMULATION OF GRAVITY AND STOCHASTIC 

DIFFERENTIAL EQUATIONS 

We shall adopt the standard Lagrangian density in Einstein's gravitational theory: 

C^Cg + Cm, (72) 

17 



where 



I^G = V^g^^'R,., (73) 

Cm is the Lagrangian density for the matter field, g = det g^u and -R^,y is the Ricci curvature 
tensor. 

The control parameter is defined in terms of spacetime fiuctuations of the gravitational 
constant: 

G, = G + aC,, (74) 

where, as in Eq. (H), G is the average value of the gravitational constant, and denotes 
Gaussian spacetime noise. The diffeomorphism invariant Einstein field equations including 
stochastic spacetime fiuctuations take the form: 

G^, = 87r(GT^, + ae.T^,), (75) 

where G^,^ is the Einstein tensor: G^j, = i?^,^ — \g^uR and R is the scalar curvature. 

The first step in formulating the Hamiltonian approach to gravity is to introduce a 
foliation of spacetime which defines a 3-geometry ^^^Q. We choose a time function t and a 
vector field such that the surfaces, E, of constant t are spacelike Cauchy surfaces with 
t'^'Vfj,t = 1, where denotes the covariant differentiation with respect to the metric tensor 
gfj^u- In contrast to Minkowski spacetime, there is no preferred coordinate system in curved 
spacetime. In the standard (3 + 1) treatment of spacetime |jl3|-|l6|, we introduce the lapse 
function, A^, by 

N = -g^.t^n" = (n^'V^t)-' (76) 

and the shift vector A^^ by 

N'' = h'^.t", (77) 
where is the unit normal vector to E, 
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(78) 



is the induced spatial metric on E, and we use the metric signature: (—1, +1, +1, +1). We 
have 



In the following, the volume element is 



The scalar curvature, i?, is given by 



We also have 



(79) 



(80) 



(81) 



(82) 



where K^i, is the extrinsic curvature of E defined by 



MJ5 



(83) 



and -D^ is the derivative operator on the surface E connected with h^u. This leads to the 
expression for Cg obtained by Arnowitt, Deser, and Misner 



The Hamiltonian density associated with Cg is given by 



where vr'^'^ is the canonically conjugate momenta to h^y defined by 
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(84) 



(85) 



(86) 



In (p5D, the and A*"^ are not treated as dynamical variables but, instead, define the 
configuration space in terms of the metric, h^^. 
Variation of 



HG = J^nG, (87) 

leads to the constraint equations: 

^^'^R - h-^TX^'^'K^^ + hi-^Tx^^,Tx\ = 167rGTx±, (88a) 

D^{h-^/^Tx^"') = -MGTl, (88b) 



where Ti_l = T^uii^n^ and = h'^ aT'^^np are components of the stress energy-momentum 



tensor T^^ for matter. 



Hamilton's first order equations in time are now obtained from Hq ||T3 
5Hg 



d,h 



t't'iii' 



2/i-i/2iv(^7r^, - ^ Vvr^) + 2D^^N^), (89a) 



5V ^ 2 J 2 V 2 

+ /i^/2^^(/^-i/2^<T^M^ _ 27r'^(/^D^A^'')) + SvrGT^^ (89b) 
We shall now write the stochastic differential equation for the dynamical random variable 



Ht as 



dtT^I" = + SnGT^" + 87^a6T'^^ (90) 

where 



1 

2 M a 



+ /i^/2D,(/i-i/2jya^A.z. _ 27r"(^D^iV'^)). (91) 

Exploiting the fact that the stochastic equation is the derivative of the Wiener process, and 
using the control parameter, Gt = G + a^ti we obtain 
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c/vrr = Fr {7Tt)dt + BrdWt , (92) 

where 

i^^W = /^(v^^) + 87^GT'^^ (93) 

and 

Br = StktT'"'. (94) 

The constraint equations ( p8a|) and (^8b|) take the stochastic form: 

- + Ik'^^K. = l67i{GT^^ + a^tT^^), (95a) 

D,ih;'/\r) = -^GTl + aitTl). (95b) 

Consider the random process 

Zt = v{Tiut). (96) 



According to the Ito rule |T^, the SDE of the process Zt is 

dZt = dv{7Ct, t) = [dtv{7c, t) + Ft{7Ct)d^vi7f, t) + ^B^d^^-^, t)]dt + Btd^v{Ti, t)dWt, (97) 

where for simphcity we have written vr^'^, F^*^, B^^'^ and T^^ as tt, F, B and T. 
The integral form of (^) is given by 

ft' 

v{7rt>,t') - v{no, 0) = / [dtv{nt, t) + Ft{nt)d^v{nt, t) 

Jo 

1 /■*' 

+ -B^d„M'^t,t)]dt+ Btd^v{Tit,t)dWt. (98) 

2 Jo 

We assume that v{n,t) has compact support, so that 

1. v{n,0) = 0, 

2. 11(71, 00) = 0, 

3. j;^ E{{Btd^v{TT,t)f }dt < 00. 
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We have 



BtdM^t, t)dWt'^ = 0, (99) 

and for t' ^ oo in (p8| ) and taking the expectation value, we get 

/ E{[dtv{7Ct,t) + Ft{-Kt)d^v{-Kut) + -Bld^^viTiumdt = 0. (100) 
Jo 2 

If the transition probabihty of the stochastic process vTf possesses a density, then we have 



dt / dypiy,t\Trt,0)[dtviy,t) + F{y)dyv{y,t) + 327T'a'T'dyyv{y,t)] = 0. (101) 
Jn 

By performing a partial integration, we obtain 

/•oo r 

/ dt / d?/t;(y,t)[-ai;?(2/,t|7r,,0)-9jF,(y)p(y,t|7ri,0)] + 327rVT2a^,p(?/,t|7r,,0)] = 0. 

(102) 

Since v{'Kt.,t) is an arbitrary function, it follows from ( |102| ) that 

dtp{y, tlvTi, 0) = -dy[Ft{y)p{y, t|7ri, 0)] + ?,2Ti^a^T^dyyp{y, t\'Kt, 0). (103) 

Thus, the probability density for the random conjugate momentum variable vrf'^ satisfies a 
Fokker-Planck equation. 

The stationary probability behavior of a gravitational system is given by 

c fir F(u) \ 

where C is a normalization constant and we have used the Ito interpretation of (pO|). 

VII. STOCHASTIC EQUATIONS OF MOTION 

Consider the geodesic equation: 

du^" 



ds ^'"^ 
where 
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+ Kpu'^u" = 0, (105) 



= \9^'' {do^dap + dpgaa ' d.Qap) , (106) 

denotes the Christoffel symbol, it^ = dx^^/ds denotes the time-hke four- velocity and ds is 
the increment of proper time along the world line, defined by 

ds^ = -g^^dx^dxr (107) 

We can then consider to be a random variable, m^, and form the SDE: 

(i< + r^,,^«>fcis = -CsF^ds, (108) 

where is a vector quantity random variable, formed from the metric tensor and its 
derivatives, is a Brownian motion process in terms of the proper time s, and the Christoffel 
symbol is treated as a random variable determined by the stochastic metric gs,fiv At the 
length scale for which the fluctuations of spacetime are significant, we picture a test particle 
moving in spacetime along a Brownian motion path such that u'^ does not have a well-defined 
derivative with respect to s at a point along the world line. For larger macroscopic length 
scales for which the spacetime fluctuations can be neglected, the motion of the test particle 
becomes the same as the deterministic geodesic motion in GR. 
Let us define the spatial metric: 

h^j,^ = g^j,^ + u^u^. (109) 

Then, hi^j, — g^h^i, is the projection operator onto the subspace of the tangent space 
perpendicular to u^. The expansion 6', shear a^i, and twist uj^^ of a congruence of geodesies 
are defined by 

= V> (110a) 
= y{f,y) - ^Oh^^, (110b) 
^^^lv^y[nv], (110c) 

where the tensor field, F^,^, is 
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;iii) 



and is purely spatial: 



(112) 



The vector field, u^, of tangents is normalized to unit length, u^u^ = —1, and is defined 
by 



(113) 



The Raychaudhuri equation takes the form |26,29,15 



flO 1 



ds 



Using Einstein's field equation, we have 



We assume that the weak energy condition: 



(114) 



(115) 



T^uU^u" > 



(116) 



holds for matter for all timelike m^. Moroever, we also assume that the strong energy 
condition holds: 



T^uU^u'' > --T\. 



;ii7) 



We choose the geodesic congruence to be hypersurface orthogonal, so that u^j^y = 0, 



whereby the third term in (|114|) vanishes. It follows from the condition ( |117|) that the left 
hand-side of ( p.l4| ) is negative. Then we have 

de 1 



(118) 



which gives 
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and, therefore, 



e-\s) > 9,' + ^s, (120) 



where 60 is the initial value of 6. If we assume that 60 is negative for a converging congruence 
of geodesies, then ( |12CI| ) implies that 6^^ must pass through zero and 6 — > —00 within a 
proper time s < 3/16*01. Thus, it follows for a tube of matter that at a point where 6 —>■ —00, 
the matter density p ^ 00 and there is a singularity at that point on the world line. 
For null geodesies, the Raychaudhuri equation takes the form: 

^ = -^e' - a,,a^^ + - R,.k^k\ (121) 

where we now consider a congruence of null geodesies with the tangent field k^. The energy 
condition takes the form: 

T^^Pr > 0. (122) 



Let us write ( |114| ) in the form 



df) 1 

_ = -(-e' + a,,a^^ + SttGT), (123) 

where we have again chosen a congruence which is hypersurface orthogonal, uj^,^ = 0, and 

f = T^,n^M^ + (124) 

We can now use the stochastic control parameter, Gs = G + a^g, and obtain the SDE: 

dOs = -{1:91 + a,,^,crr + SnGf)ds - SnafdWs, (125) 

where dWg denotes the Wiener process, dWg = dC,sds. For null geodesies, we obtain the 
SDE: 

dOx = -{^el + (TA.^^af + 87rGf)dX - SnafdWx. (126) 
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We see that for a given length scale if the intensity of the fluctuations is big enough, 
then the left hand-sides of ( |125| ) and ( |126| ) are no longer negative definite and it no longer 
follows that 1^1 ^ oo at points along the world line. Caustic singularities can be prevented 
from developing in a congruence of timelike or null geodesies if convergence occurs anywhere 
in the manifold, provided the geometrical fluctuations are strong enough. Thus, stochastic 
gravity can avoid the occurrence of singularities in spacetime. In the limit of classical GR, 
the Hawking-Penrose theorems hold, for the Brownian motion fluctuations of spacetime are 
negligible and can be neglected. 

The big-bang singularity in cosmology and the singularity which occurs in the gravita- 
tional collapse of a star, in GR, can both be avoided in stochastic gravity. It is to be expected 
that in the limit of zero spatial dimensions and zero time, the spacetime fluctuations will 
dominate and smear out singularities. In the next section, we shall see that it is also possible 
to avoid the occurrence of infinite red shift event horizons in the collapse of stars. 



VIII. GRAVITATIONAL COLLAPSE 

Let us apply our stochastic gravitational theory to the problem of gravitational collapse 
of a star. A cooling star of mass greater than the Chandrasekhar or Oppenheimer-Volkoff 
mass limit cannot maintain equilibrium as either a white dwarf of a neutron star, and is 
predicted in GR to collapse to a black hole. We shall restrict ourselves to the case of 
spherically symmetric collapse of "dust" with negligible pressure. In practise, we should 
solve Hamilton's first order equations ( |89a|) and (|89b| ) , since ( ^9b|) is linear in the control 



parameter G and we can solve the stochastic equation ( p2D directly in the Gaussian white- 
noise limit. However, instead we shall solve the collapse problem by following the treatment 
given by Oppenheimer and Snyder and by Weinberg P7|J25[] . 

We use a comoving coordinate system to describe the freely falling dust particles. The 
metric is given in Gaussian normal form by 

ds^ = - U{r,t)dr^ - V{r,t){de^ + smHd(j)^), (127) 
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while the energy-momentum tensor for the fluid is 



;i28) 



where p(r, t) is the proper energy density and is given in comoving coordinates by 



= u'' = uf = 0, = 1. 



(129) 



The energy conservation equation is 



dtipVVu) = 0. 



(130) 



The Einstein field equations are given by 



131) 



where 



(132) 



Assuming that the collapse is homogeneous, we can seek a separable solution: 



U = R\t)f{T), V = S\t)g{T). 



(133) 



Then, Einstein's field equations require that |^ S/S = R/R, where R = dtR. We can 
choose : S{t) = R{t) and redefine the radial coordinate, so that V = R^{t)r'^. Solving the 
resulting Einstein field equations leads to the Friedmann- Robert son- Walker (FRW) metric: 



ds' = dt' -R\t) 



.1 — kr"^ 

We normalize R{t) so that -R(O) = 1 and obtain: 



+ r'^de'^ + sin^ 



(134) 



p{t)=p{m-\t). 



(135) 



The field equations yield 
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(136) 



where A; is a constant given by 



k = ^p(O). (137) 



Eg. ( |136| ) becomes 

1 



R' = k 

with the parametric cycloid solution: 



R{t) 



1 



;i38) 



i? = i(l + cosV^). (139b) 



At time t = ts, where 



the fluid sphere collapses to an infinite proper energy density. The resulting singularity as 
the end stage of collapse is inevitable in classical GR, as follows from the Hawking-Penrose 
theorems on gravitational collapse ||29 



According to the Birkhoff theorem, the metric outside the collapsing star is static and 
given by the Schwarschild solution: 



2GM\ , , / 2GM\-i , , 



1-— — dr^ -r\de^ + smHd(P^). (141) 



By choosing an integrating factor, we can transform the metric to a standard form p8| : 

ds'^ = B{r, t)dt^ - A(r, t)dr'^ - r^{de^ + sin^ Odcf)'^), (142) 

where 

^-^It^J (1 - kryR) ' ^^^^^^ 
^=(l--) . (143b) 
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and the constant a is equated to the radius of the sphere in comoving polar coordinates. It 
follows that the interior and the exterior solutions match continuously at r = aR{t) when 



k ^ ^, (144, 



a3 



which gives M = ^p{0)a^. 
Let us rewrite Eq. (|138|) as 



/2GM\i/2 



1 



1 



1/2 



(145) 



where for the collapse problem we have chosen the negative square root. The nonlinear SDE 



has the form: 



dRt = [VGfiRt) + P{VG + Ct)f{Rt)]dt, (146) 



where 



For rapidly varying geometrical fluctuations, corresponding to a random process with a C,t 
that has a short correlation time, we can obtain in the white-noise limit the approximate 
stationary probability density: 



[R dR 



J v{R) 



148) 



where C is a normalization constant and 

v{R) = f{R)[l + T,MR) + ...]. (149) 

To lowest order we get 

We see that i? = 1 is a natural boundary [|T7|: both the drift and the diffusion coefficients 
vanish at i? = 1. For \/G < (T^/2(a^/2M)^/^ the point i? = 1 is attracting and the 
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stationary probability is mostly concentrated at i? = 1. The point a/G = cr^/2(a^/2M)^/^ 
is a transition point, since for > aV2(aV2M)i/2 ^j^g probability density moves away 
from the divergent point R = 1 and a new stationary probability density emerges with a 
large finite probability for collapse towards R = 0. However, for R ^ 0, we have 

PsiR) ~ 0. (151) 

We therefore arrive at the result that as the star collapses there is zero probability for R{t) 
to have the value zero, and consequently there is zero probability of having a singularity as 
the final state of collapse. 

Let us now consider the red shift emitted from the surface of the star as it collapses. The 
fractional change of wave length emitted at the surface is 

where 

t' = tU' (l-?£^y,r (153) 

JaR{t) \ r J 

is the time it takes for a light signal emitted in a radial direction at standard time t to reach 
a distant point r. In the limit, R{t) IGMja = ka"^, we obtain by using (|145|) : 

We have 

dz = f{G,R)dt, (155) 



where 



We can now define the nonlinear SDE: 
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dzt = ^i{G)u{Rt)dt + 7(G + Ct)u{Rt)dt, (157) 

where the function u{R) oo as R ^ 2GM/a, corresponding to the infinite red shift hmit 
as the event horizon is approached during collapse predicted by classical GR, and Q is a 
short time stochastic, colored process. Expanding around the white-noise Gaussian limit, 
we can derive the approximate stationary probability density for the red shift in the limit 
R 2GM/a: 



C 

u^{R) 



2/i(G) f dR 



J u{R). 

The result emerges from ( |158| ) that provided the bound 



(158) 



dR cr^ , / N / s 

is satisfied, then 

ps{z) ~ (160) 

and the probability for an infinite red shift at i? = IGMja is zero. Thus, in our stochastic 
gravitational theory the probability for a black hole event horizon with an infinite red shift 
to form during gravitational collapse is zero when ( |159| ) is satisfied. On the other hand, 
once a trapped surface forms in classical GR collapse, then the red shift seen by a distant 
observer must inevitably become infinite at the radius tq = 2GM. We recall that for 
physical collapse in nature, the unphysical regions which occur in the Penrose diagram for 
the extended (Kruskal) Schwarzschild solution are covered up by the collapsing matter. In 
GR, the metric outside the collapsing star must be the Schwarzschild static metric (Birkhoff 
theorem) . 

In stochastic gravity, the geometrical fiuctuations about the deterministic event horizon 
cut off the high wave length at a finite value, A = Ac, when ( |159| ) holds, as viewed by an 
observer at large distances from the collapsing star. In fact, the distant observer does not see 
a collapsed object described by a static Schwarzschild metric with no hair. The collapsed 
object will fluctuate about the average Schwarzschild solution with a characteristic time 
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determined by the size of the correlation function for the Schwarzschild metric. The red 
shift of the collapsed object can be high, so that an outside observer would believe that 
it is a black hole. The cooperative effects associated with the self-organising microscopic 
subsystems comprising the event horizon can produce a large or even infinite correlation 
length for the microscopic fluctuations, cutting off the infinite wavelength radiation emitted 
by the macroscopic surface of the collapsing star as R ^ 2GM/a. 

Our use of the stationary probability density for the collapse problem means that we have 
assumed that after an infinite time has elapsed the system evolves to a stationary state. This 
is compatible with the fact that the spherically symmetric collapse problem has well-defined 
fiat space asymptotic limits as r — oo. Since in stochastic gravity, it can take a long (but 
finite) time for a distant observer to see the formation of an apparent event horizon during 
collapse, this is consistent with the observer measuring the stationary probability density 
Ps as the limit p{z) — > ps{z) as t — oo. 

We have succeeded in removing both the singularity in collapse and the black hole event 
horizon in our stochastic probability gravitational theory (provided ( |159| ) is satisfied.) The 
deterministic GR predictions are clearly obtained in the limit that the intensity of fluc- 
tuations vanishes. As we have seen in the last section, the results are consistent with the 
Hawking-Penrose theorems for collapse, because for large enough fluctuations of the stochas- 
tic metric, gs,^ij, caustics can be prevented from occurring for converging congruences of 
timelike or null geodesies in the spacetime manifold. 

The 'greyness' of the dense collapsed object depends on the size of the correlation time 
Tcor- The correlation length can be significantly larger than the Planck length scale, Ip, since 
the fluctuations are associated with a stochastic property of the gravitational field, which is 
not directly proportional to Planck's constant h, and is therefore not strongly correlated to 
regions at the scale of Ip. If the dense collapsed object is sufficiently grey, then this could 
solve the information loss problem first introduced by Hawking Again, this depends 
on the size of the correlation length associated with the geometrical fiuctuations. Since the 
deterministic event horizon has been removed, either there is no Hawking radiation emitted 
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from the surface of the collapsed object, or the spectrum of the radiation is significantly 
modified from its behavor in deterministic GR. This raises the question: will black holes form 
in the collapse of astrophysical objects in our statistical mechanical treatment of gravity? If 
the answer is no, then what kind of object would describe the final state of collapse? These 
interesting questions require further investigation, before any conclusions can be reached 
about this fundamental problem. 

IX. THE COSMOLOGICAL BIG-BANG SINGULARITY 

Let us now apply our stochastic formalism to cosmology. We shall consider the standard 
FRW cosmology. The FRW metric has the form ( |134| ) and solving Einstein's field equations 



including the cosmological constant term yields the Friedmann equation: 

where R{t) denotes the cosmic scale factor, k now denotes the measure of the spatial cur- 
vature: k = 1,0,-1, po^M and po.r denote the present density of matter and radiation, 
respectively, and A denotes the cosmological constant. 

Let us consider for simplicity the case of zero spatial curvature and zero cosmological 
constant, A = A; = 0. We shall focus on the radiation dominated universe near t = 0. Then, 



Eq.( |161j ) can be written as 

ifl=yG(5^)'"ijdi. (162) 

We can now form the nonlinear SDK: 

dRt = VGf{Rt)dt + a{VG + Ct)f{Rt)dt, (163) 

where Q is assumed to be a short correlation time stochastic process associated with colored 
fluctuations of the geometry, and 



m) = ['-^) ^. (164) 



1/2 I 
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The approximate stationary probability density calculated in the limit of white-noise is 
given by 



ps{R) ~ C[ ) exp 



(165) 



cr^ / V247rpo,r- 

where C is a normalization constant. We have i? — as t ^ and we see that the 
probability of R reaching t = is zero. Thus, in our stochastic gravity the probability 
of the big-bang singularity occurring is zero. As in the case of gravitational collapse, the 
fluctuations of the geometry near t = smear out the singularity. The result is consistent 
with the Hawking-Penrose theorem [^], which states that in classical GR, the universe 
must collapse to a singularity at t = provided the positive energy conditions are satisfied, 
because for large enough fluctuations in the neighborhood of t = 0, the converging geodesic 
congruences can be prevented from forming a singularity, as has been demonstrated in section 
VII. 



X. CONCLUSIONS 

We have formulated a self-consistent gravitational theory based on stochastic diffusion 
processes. The geometry of spacetime is assumed to be like a fluctuating environment with 
the control parameter, = G + a^^. The fluctuations of the control parameter form a 
bell-shaped distribution about the deterministic Newtonian gravitational constant G. This 
theory represents a statistical mechanical treatment of gravity with a clearly determined 
classical GR limit when the intensity of the spacetime fluctuations vanishes. In this ap- 
proach, Einstein's gravitational theory is a deterministic "effective" theory, which holds to a 
very good approximation for macroscopic gravitational systems with scales larger than the 
correlation lengths associated with the geometrical fluctuations. The stochastic gravity the- 
ory does not represent a complete quantum mechanical theory and cannot solve the problem 
of quantum gravity, but it does lead to a methodology that can account for the complex- 
ity of gravitational phenomena as are experienced in cosmology, gravitational collapse and 
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astrophysics. Since it is a probabilistic theory, it reconciles a probabilistic interpretation 
of the energy-stress tensor of matter and the geometry of spacetime. It leads to critical 
phenomena and self-organization of gravitational systems, which are strongly expected to 



play a significant role in astrophysics and cosmology |31[ . 

By working with a (3 + 1) foliation of spacetime, and using the Arnowitt, Deser and 
Misner formalism, we were able to obtain a stochastic differential equation for the random 
canonically conjugate momentum vrf^. We also derived a Fokker-Planck equation for the 
probability density of vTf"^ and showed how a stationary probability density for the dynamical 
variable tt^'^ could be obtained. 

Stochastic geodesic congruences were studied and from the Raychaudhuri equation, we 
derived a stochastic differential equation for the expansion parameter 6, which led to the 
result that, assuming the weak and strong positive energy conditions are satisfied and for suf- 
ficiently strong spacetime fluctuations at small scales, a congruence of converging geodesies 
can be prevented from forming caustic singularities in the manifold, thereby avoiding the 
classical GR Hawking-Penrose singularity theorems. 

For certain applications of the formalism, we found it necessary to consider nonlinear 
non-Markovian colored processes which were expanded about the white-noise Gaussian limit, 
a limit expected to hold for a rapidly fluctuating spacetime geometry. Stationary state 
probability densities were obtained in the Gaussian limit with a short correlation time for 
the fluctuating geometry in gravitational collapse and cosmology. It was found that the 
fluctuations of the geometry smeared out the singularity in gravitational collapse. Similarly, 
the big-bang singularity was smeared out in a probabilistic sense, leading to a singularity-free 
FRW universe. 

A stochastic treatment of gravitational collapse resulted in the removal of black hole 



event horizons (given that the bound ( |159| ) is satisfied) by cutting off the high frequencies 
in the infinite red shift limit. This means that either there is no Hawking radiation emitted 
from the surfaces of the collapsed objects, or the blackbody radiation spectrum predicted 
by deterministic GR is significantly modified. Moreover, the lack of a rigorous event horizon 
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could remove the problem of information loss [0 depending on the 'greyness' of the collapsed 
objects. 

We can also extend the stochastic methods developed here to the nonsymmetric grav- 
itational theory (NGT) |]32[, so that for small enough fluctuations of the non-Riemannian 
geometry we obtain the classical NGT theory. 
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